را فرش 3 ات :01 سال سی و دوم: شمارٌ یک» ۱۳۹۸ 


مطالعة عددی اثر زاوية اریب در کانال م رکب با سیلاب‌دشت‌های مورب بر اند ر کنش بین جریان در کانال 
اصلی و سیلاب‌دشت‌ها* 
(یادداشت پژوهشی) 


محمدمهدی سیف( بهرام رضائی( 

چکیده در این تحفیق تلاش شده‌است که میدان جریان شامل توزیع سرعت و تنش برشی مرزی د رکانال‌های مرکب با سیلاب‌دشت‌های 
مورب برای دو زاویهٌ اریب ۵/۱ و ۹/۲ درجه و عمق‌های نسبی ۰/۱۵ ۰/۲۶ ۰/۶۱ و ۰/۵۰ به صورت عددی شبیه سازی شود. به این منظور از 
مد ل 7 شفتکی ۲-۶ و نرم‌افزار ۱۸۷۹۳-0۳۲ ستفاده شدها ست. ننایج حا صل از شبیه سازی عددی توزیع سرعت و تتش بر شی در دو مقطع 
انتحاب و با داده‌های به‌دستآمده د رآزمایشگاه هیدرولیک مسقر در والینگفورد (وذان 60۵ ۳۱000) کشور انگلستان مقایسه شده‌است. 
بررسی‌ها نشان‌دهنده انطباق قابل قبول بین داده‌ها ی آزمایشگاهی با نتایج حاصل از شبیه‌سازی عددی می‌باشد. با کمک معادلات اندازه حرکت 
توسعه داده‌شاه برای حجم کنترل‌هایی در سیلاب‌دشت‌هاء اند رکنش بین جریان د رکانال اصلی و سیلاب‌دشت‌ها برای زوایای اریب و 
عمق‌های نسبی مختاف, مورد مطالعه قرا رگرفته‌است. بررسی‌ها نشان‌دهند ‏ آن است که عموماً به‌دلیل تبادل جرم بین کانال اصلی و 
سیلاب‌د شت‌ها با افزايش زاویهة اریب و عم قآب, نیروی بر شی ظاهری در فصل مشترک قائم بین جریان د رکانال اصلی و سیلاب‌د شت‌ها 


افزایش می‌یابد. 


واژه‌های کلیدی شبیه‌سازی عددی. کانال مرکب با سیلاب‌دشت‌های مورب. مدل آشفتگی - نیروی برشی ظاهری. 
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تقعطه اجم‌ندم۸ راممجظ عمصماناطایتط عی! روامصصقطه مصتممرومی 9۵0ما رعطتاممصر لهمتمصت ۲۷۷۵۲۵5 رمک 
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* تاریخ دریافت مقاله ٩0/۷/۲۳‏ و تاریخ پذیرش آن۹۸/۳/۲۹ می‌باشد. 
(۱) فارغ‌التحصیل کارشناسی ارشد عمران» مهندسی آب و سازه‌های هیدرولیکی, دانشگاه بوعلی سینا. 


(۲) نویسنده مسئول استادیار گروه عمران دانشکده مهندسی, دانشگاه بوعلی سینا. 028۱26 ۵) [۲۵226 ۲۱۵2۲1:۵۰ 


۱۲ 


مقدمه 
به‌ندرت یافت می‌شود. شکل مقطع و پروفیل طولی 
رسوبت در طی قرن‌ها شکل داده ستله آنینتا: بنابراین 
رودخانه‌ها دارای شکل مقطع عرضی کاملاً پیچیده با 
ها همواره به دلایل مختلف نظیر داشتن خاک حاصل - 
خیز و نزدیکی به ساحل رودخانه مورد توجه بشر قرار 
داشته و بسیاری از تمدن‌ها در ساحل رودخانه‌های 
بزرگ شکل گرفته‌اند؛ بنابراین آشنایی با هیدرولیک 
جریان در سیلاب‌دشت‌ها برای تأمین حفاظت جان 
انسان‌ها و نیز ساژه‌ها و تأسیسات موجود فر این 

نواحی لازم و ضروری است. 
در کانال‌های مرکب. توزیع عرضی سرعت تحت 
تأثیر بالا بودن ضریب زبری بستر و عمق کم جریان 
به‌شکل قابل ملاحظه‌ای حالت 
غیریکنواخت دارد. به‌دلیل اختلاف سرعت بین جریان 
در کانال اصلی و سیلاب‌دشت ها؛ یک لاه برشی در 


سیلاب دشتی» 


فصل مشترک این نواحی تشکیل می‌شود و در نتیجة 
آن گردابه‌های بزرگ‌مقیاس با محور قائم و نیز جریان 
ثانویةٌ مارییچی شکل درجهت طولی ایجاد می‌شود 
(تومیناگا و نزو [1]). این در حالی است که در 
کانال‌های مرکب با سیلاب‌دشت‌های مورب به‌دلیل 
تبادل جرم و اندازهٌ حرکت از سیلاب‌دشت‌ها به کانال 
اصلی و برعکس. آشفتگی‌های با مقیاس به‌مراتب 
بزرگ‌تر در فصل مشترک کانال اصلی و سیلاب‌دشت‌ها 
به‌وجود می‌آید. 

جیمز و بارونز [2] نخستین محققانی بودند که به 
بررسی آزمایشگاهی جریان در یک کانال مرکب با 
کانال اصلی مورب نسبت‌به سیلاب‌دشت‌ها پرداختند. 
آنها دریافتند که به‌دلیل تبادل جرم. شتاب جریان در 
سیلاب‌دشت واگرا افزايش می‌یابد درحالی‌که در 


نشريه مهندسی عمران فردوسی 


مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 


سیلاب‌دشت همگرا کاهش می‌یابد. 

الیوت [3]. الیوت و سلین [4] و سلین [5] 
بااستفاده از تجهیزات مستقر در والینگفورد ۳0۳ 
(رانانه‌۲ اعمعدط) ۳۱۵۵۵) ۰ به‌صورت آزمایشگاهی 
جریان در کانال‌های مرکب با سیلاب‌دشت‌های مورب 
را مورد بررسی قرار دادند. آنها دریافتند که توزیع 
سرعت و تنش برشی مرزی در کانال اصلی و سیلاب- 
دشت‌ها به‌شدت تحت ۳3 جریان عرضی مبادله‌شده 
بین این نواحی قرار دارد. آنها هم‌چنین بااستفاده از 
معادلات انداز حرکت. نیروهای عمل‌کننده بین کانال 
اصلی و سیلاب‌دشت‌ها را مورد مطالعه قرار دادند. 
چلبک [6] داده‌های به‌دست‌آمده توسط الیوت [4] و 
سلین [5] را مجدداً تجزیه‌وتحلیل نمود و با نتایج 
تست آ لو ۱۰ فلوم میرن متفر قن ذانشگاه 
بیرمنگهام مورد مقایسه قرار داد. 

در رابطه با شبیه‌سازی عددی پزینگا [7] بااستفاده 
از مدل غیرخطی ع- به شبیه‌سازی میدان جریان در 
یک کانال مرکب منشوری پرداخت. او دریافت که مدل 
آشفتگی مورد نظر قادر به پیش‌بینی جریان ثانويه ناشی 
از تتش آشفتگی غیرایزوتروپ است. بسیاری از 
محققان مانند وو و همکاران [8]» فیسچر- آنتز و 
همکاران [9] و راماش‌واران و نادن [10] برای 
شبیه‌سازی سه‌بعدی جریان در کانال‌های باز مدل 
آ مشک ۶-۵ را ترجیح داده‌اند. کانگ و چوی [11] به 
مطالعهٌ جریان در یک کانال مرکب منشوری بااستفاده 
از مدل آشفتگی تنش رینولدز 50699 ۴680109) 
([1000 پرداختند. آنها برای راستی‌آزمایی نتایج 
شبیه‌سازی عددی, از داده‌های تومیناگا و نزو [1] 
استفاده کردند. 

کوکالجات [12] جریان در کانال‌ها را بااستفاده از 
دو مدل آشفتگی ع- و مدل تنش انتقالی رینولدز 
(اع0مصه )128۵0۲ ووعع9 ۵78۵108) شبیه‌سازی 
کرد. او دریافت که مدل تنش انتقالی رینولدز قادر است 


که موقعیت و قدرت سلول‌های جریان ثانویه و نیز 


سال سی و دوم شماره یک» ۱۳۹۸ 


محمدمهدی سیف- بهرام رضائی 


نقش آن را در جابه‌جایی موقعیت سرعت حداکثر به 
ناحیه‌ای در زیر سطح آب» پیش‌بینی کند. اما در مقابل 
مدل‌های غیر حطی 8 قادر به انجام اين امر نمی - 
بان 

رایت و همکاران [13] مدل آشفتگی ع- و 
مدل‌های مختلف تنش رینولدی را برای مطالعة 
مشخصات جریان در یک کانال با مقطع عرضی 
ذوزنقه‌ای به کار برد. نتایج این بررسی نشان داد 
درحالی که همه مدل‌های آشفتگی مشخصات کلی 
جریان را به‌شکل یکسان و مشابه پیش‌بینی می‌کنند. ما 
نتایج آنها در شیبه‌سازی جریان انویه کاملاً متفاوت 
است. 

بیمن [14] بااستفاده از مدل آشفتگی شبیه‌سازی 
گردابهٌ بزرگ (طمنامابهنگ 50۵0 6عتض) به شبیه- 
سازی جریان در کانال اصلی و جریان سیلاب‌دشتی 
پرداخت. نتایج شبیه‌سازی گردابة بزرگ نشان داد که 
این مدل با دقت بالایی قادر به شبیه‌سازی جریان‌های 
انویه در جریان‌های درون کانال اصلی و خارج آن 
(جریان سیلاب‌دشتی) می‌باشد. 

رضائی و امیری [15] به شبیه‌سازی جریان در 
کانال مرکب با سیلاب‌دشت‌های همگرا پرداختند. 
بدین منظور آنها با کمک نرم‌افزار ۸۲۹۷5 و بااستفاده 
از مدل آشتفگی عع توزیم سرعت و تنش برشی را 
برای کانال مرکب غیرمنشوری با دو هندسه متفاوت 
مورد بررسی و با داده‌های آزمایشگاهی مورد مقایسه 
قرار دادند. نتایج نشان داد که مدل آشفتگی 1-6 
توانایی شبیه‌سازی توزیع سرعت و تنش برشی را با 
دقت قابل قبولی دارد و نیز می‌تواند وجود سلول‌های 
جریان انویه در کانال اصلی ناشی از انتقال جرم از 
سیلاب‌دشت‌ها را به خوبی پیش‌بینی نماید. 

در این تحقیق سعی بر این شدها ست که میدان 
جریان در کانال مرکب با سیلاب‌دشست‌های مورب 
(زوایای اریب ۵/۱ درجه و ۹/۲ درجه) به کمک مدل 


آشفتگیع-* شبیه‌سازی و اندرکنش جریان بین کانال 


مان سی و دوم شماره یک. ۱۳۹۸ 


۱5۳ 


اصلی و سیلاب‌دشت‌ها مورد مطالعه قرار گیرد. 
معرفی کانال مرکب با سیلاب‌دشت‌های مورب 
کانال مدل‌شده مشابه کانال آزمایشگاهی ۳0۳ مستقر 
در والینگفورد و مورد استفاده به‌وسیلهٌ الیوت و سلین 
[4] است که دارای شیب ۰۱/۰۲۷*۱۰۲ طول ۵۵ مترء 
عرض ۹/۰ متر و از جنس بتنی صاف است. به‌منظور 
انجام آزمایش برروی کانال مرکب با سیلاب‌دشت‌های 
مورب مقطع عرضی به‌شکل مرکب ذوزنقه‌ای‌شکل با 
کانال اصلی به‌عرض ۱/۵ متر. عمق ۰/۱۵ متر وشیب 
دیوارهٌ ۱:۱ تغییر شکل داده شده‌است. هم چنین برای 
ایجاد سیلاب‌دشت‌های مورب. دیواره‌های متحرک در 
کانال نصب شده‌است. شکل‌های (۱) و (۲) تصویر 
عمومی. نمای شماتیک و مقطع عرضی کانال 

آزمایشگاهی مورد استفاده را نشان می‌دهند. 

اندازه گیری سرعت نقطه‌ای در دو مقطع شماره ۱ 
و ۲ به فواصل ۱۰ متری از یکدیگر بااستفاده از 
سرعت‌سنج پروانه‌ای با قطر خارجی ۱۶ میلی‌متر» در 
فواصل عرضی به‌ترتیب ۵۰ و ۱۰۰ میلی‌متری در کانال 
اصلی و سیلاب‌دشت‌ها و نیز فواصل ارتفاعی ۱۰ 
میلی‌متری انجام شده‌است. هم‌چنین به‌منظور برآورد 
اندازژ حرکت مبادله‌شده بین کانال اصلی و سیلاب- 
دشت‌هاء در فواصل طولی ۱ متری. مقادیر سرعت در 
فصل مشترک بین کانال اصلی و سیلاب‌دشت‌ها 
اندازه گیری شده‌است. 

تنش برشی موضعی مرزی نیز بااستفاده از لولة 
پریستون به‌قطر خارجی ۶/۰۲ میلی‌متر در پیرامون 
محیط تر شده کانال با فاصله‌های ۵۰ میلی‌متری برروی 
دیوارها و نواحی نزدیک دیواره‌ها و ۱۰۰ میلی‌متری در 
کف کانال اصلی و سیلاب‌دشت‌ها اندازه‌گیری 
شده‌است. این اندازه گیری در همان مقطعی که سرعت 


متوسط در عمق به‌دست مک هه انجام شده است. 


نشریه مهندسی عمران فردوسی 


۱5 


شکل ۱ تصویری از کانال مرکب با سیلاب‌دشت‌های مورب 


(زاوية اریب ۵/۱ درجه) مستقر در والینگفور 


صنجام۲[000 
68665 0۲9.2 5.1 


ا<,0 


۵1 


شکل ۲ نمای بالا و مقطع عرضی کانال مرکب با 


سیلاب‌دشت‌های مورب 


دینامیک سیال محاسباتی یک ابزار محاسباتی 
قوی برای شبیه‌سازی میدان جریان در سیالات و پدیده 
انتقال حرارت در فرایندهای فیزیکی است. این روش 
با حل عددی معادلات متوسط گیری‌شدهة ناویر- 
استوکس در یک میدان مشخص با شرایط مرزی 
توت شاه عمل ع ی کل 

میدان جریان مدل‌شده در این تحقیق. کانالی 
است به‌طول ۵۵ متر و عرض ۹/۱ متر. عرض کانال 
اصلی ۱/۵ متر و سیلاب‌دشت‌های مورب با دو زاويه 
اریب ۵/۱ و ۹/۲ درجه. شبیه‌سازی برای دبی‌های 


نشريه مهندسی عمران فردوسی 


مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 


۰ ۰/۳۹۰ ۰/۷۰۰ و ۱/۱۰۰ متر مکعب بر ثانیه 
متناظر با عمق‌های نسبی (01(/۴)-۳) ۰/۱۵ ۰/۲۶ 
۱ او ۰/۵۰ انجام شده و با داده‌های 
آزمایشگاهی ۳0 سری ۸14 و ۸15 مقایسه شده‌است 
(در اين رابطه 11 عمی آب و " عمق کانال اصلی 
نسبت‌به سیلاب‌دشت‌ها است). 
شبیه‌سازی عددی جریان در کانال مرکب با 
سیلاب‌دشت‌های مورب 
معادلات حاکم 

معادلات حاکم بر جریان سیال شامل معادلات 
اتکی ور آنداره هر کت اسست.. این صادلانته پر آع 
جریان أشفتة ماندگار سیال تراکم‌ناپذیر به‌صورت 
تانسوری زیر نوشته می‌شوند: 


)۱( 0 - فك 


(رتارن2 .. ._ سب ۵ . نلا02 0۳5 
)۲( و2 ۵ < )وج + 2 ۳ ۳ 


پرای در نظر گرفتن اثرات مربوط به آشفتگی 
روی مولفه‌های سرعت و فشار» جریان به‌صورت 


مجموع مولفه‌های متوسط و نوسانی جای‌گذاری 

شده‌اند. 

(۳) لا ِ لا < لا 

(( ۲4 < ۲ 
در این معادلات زل] مولفة سرعت متوسط. ۴ فشار 

نام تنش‌های رینولدز است و بااستفاده از مدل 


اشفدگن لته گرذابه‌ای. پوسیتیسکه به‌ضوزت: زیر 
محاسبه هن کرهن 


2 از 01 رح 
۵( 0166 < سس ِِ با < نام 


02*1 


در معادلهة فوق نیز لا لزجت اقشتکی» 6 ار زک 
بختیشی. اتکی و .0 دکاین. کرو نکر می‌باشد: 


سال سی و دوم» شمارة یک» ۱۳۹۸ 


محمدمهدی سیف - بهرام رضالی 


معادلات ناویر- استتو کنبن به‌صورت 
متوسطگیری‌شده در زمان یک دستگاه معادلةً 
دیفرانسیل متشکل از چهار معادله با ده مجهول (شامل 
چهار مجهول :لاو ۴ و شش مجهول تنش‌های رینولدز 
(ایدام-) است. بااستفاده از مدل آشفتگی لزجت 
گردابه‌ای بوسینیک شش معادلهٌ اضافی به‌دست می‌آید 


و دستگاه معادلات تکمیل می‌شو د. 


معرفی مدل آشفتگی ع-]1 

و این سل ۴ اترزی نختیشی آشفتگی انست و فاریاس 
نوسانات در سرعت تعریف می‌شود و ۶ پراکند گی 
(اتلاف) آشفتگی گردابه‌ای می‌باشد (نرخی که در آن 
نوسانات سرعت تلف ی کرک مدل ع-1 دو متغیر 
جدید در سیستم معادلات وارد می‌کند. در مدل 1-2 
فرض می‌شود که لزجت آشفتگی به انرژی جنبشی 
آشفتگی و اتلاف انرژی آشفتگی از طریق معادلةٌ زیر 
مربوط می‌شود: 
0( 6 از 

رتیت سیک ی بر 294 سم با ناکت 
و 6 به صورت مستقیم از معادلات دیفرانسیل انتقال 
نرخ انرژی جنبشی آشفتگی و نرخ پراکندگی آشفتگی 
بهوشت اف آزان: 


01 0 0 
۹2 ۶ ۲ + كِ 3 0 ِ (طرلام) وچ 


1 


6 ۱ 9 2 
اج( + )| - («رنام) چ 


(8موو - ببو6) ۶+ 

(۸) 

در روابط فوق ,ن» 2 6 و .6 ثابت‌هایی 
هستند که به‌ترتیب برابر ۰۱/۶۶ ۰۱/۲ ۱ و ۱/۳ 
می‌باشند. ۳۰ نیز معرف تولید آشفتگی به‌وسیلة 
نیروهای لزجت و شناوری بوده و به‌شکل زیر مدل 
ی شود: 
)۹( 


20 ازتاة ‏ ناه 
0 


9, - ۷ ) 


0 021 


سال سی و دوم» شماره یک. ۱۳۹۸ 


۱5۵ 


مشخصات شبکه‌بندی انجام‌شده 

هندسة کانال بااستفاده از نرم‌افزار 10۳102 ساخته 
و شبکه‌بندی شده‌است. برای شبکه‌بندی کانال 
وود تن از کبفاه ادن خیل ده سک ی سر 
این شده که در ایجاد شبکه المان‌های مربوط به دیوارة 
کانال اصلی و سیلاب‌دشت‌ها و هم‌چنین فصل مشترک 
بین کانال اصلی و سیلاب‌دشت‌ها به‌دلیل تغییرات 
شدید سرعت. ریزتر و در سایر نواحی درشت‌تر 
الشغات: شون 

برای بررسی تأثیر اندازةُ شبکه برروی شبیه‌سازی 
جریان سه نوع شبکه‌بندی ریز متوسط و درشت برای 
کانال مرکب با سیلاب‌دشت مورب (زاوية اریب ۵/۱ 
درجه) ایجاد شد. ابعاد عناصر سه نوع شبکه طوری 
افتسات. شدءاسته که سست انا اه ۱۳ شاد 
توسط سلیک و همکاران [16] کوچک‌تر نباشد. در 
شکل(۲) نمای شماتیک سه‌بعدی از شبکه‌بندی نشان 
داده شده‌است. بااستفاده از مدل آشفتگی 1-2 و دبی 
۷۱ متر مکعب در انیه. سرعت متوسط در عمق در 
مقطع شمارة ۱ برای هر سه نوع شبکه‌بندی درشت. 
متوسط و ریز به‌صورت عددی شبیه‌سازی و با یکدیگر 
مقایسه شده‌اند (شکل ۶). همان‌طور که در شکل دیده 
می‌شود با افزایش تعداد گره‌ها از ۲/۲ میلیون (شبکةٌ 
درشت) به تقریباً ۲ میلیون (شبکة متوسط) نتایج شبیه - 
سازی به‌صورت ملموسی تغییر نمی‌کند؛ ولی با افزایش 
تعداد گره‌ها به ۳/۹ توزیع سرعت متوسط در عمق به 
مقدار خیلی جزئی در کانال اصلی کاهش و سیلاب- 
دشت سمت راست (سیلاب‌دشت همگرا) افزایش می- 


یاید. 


شکل ۳ نمای شماتیک سه‌بعدی شبکه 


نشریه مهندسی عمران فردوسی 


۳ 


3 120 
عصز / 2 
۲-۳ 
6 ) -- 5 5ج ۱ ت‌ 
1 0 اج 
0 ص«(< 
0 ز 
(3 
0 << 
0.00 


1.4 4 


-2.6 -1.6 -0.6 4 


(حظ) 2۲ ,866 1215)2 


شکل ۶ مقایسهة سرعت متوسط در عمق شبیه‌سازی‌شده در 
مقطع شمارة ۱ برای سه نوع شبکة درشت. متوسط و ریز 


به‌منظور مقايسة نتایج شبیه‌سازی با داده‌های 
آز مایشگاهی. پراکندگی زوج مر تب های داده های 
آز مایشگاهی و عددی در اطراف خط ایده‌آل 7-6 
ترسیم و در شکل‌های (۵) نشان داده شده‌است. 
بررسسی شسکل‌ها نشسان می‌دهد که با افزایش تعداد 
گره‌ها و درنتیجه کاهش ابعاد شکه پراکندگی نقاط 
در اطراف خط ایده‌آل در ابتدا به‌مقدار جزئی کاهش 
ولی با ریزتر کردن شم که این پراک ندگی افزایش 
می‌یابد که به مفهوم کاهش دقت شبیه‌سازی عددی 
است. 

از طرف دیگر افزایش تعداد گره ها به مفهوم 
طولانی شدن زمان اجرای برنامه نیز خواهد بود. 
بنابراین به‌منظور افزايش دقت شبیه‌سازی عددی و 
کاهش زمان اجراء شبکه با انداز؛ متوسط انتخاب 


شده‌است. 


جدول ۱ جزئیات شبکه‌بندی متوسط به‌کاررفته برای 


شبیه‌سازی عددی 


ناحیه خر | تن رن ۶و9 
2 2 
۳۹ .۷ 1/۰ ۰/۱ 
8 ۲ ۱ ۱ 
‌ ۲ ۰۲۳ ۰.۲ 
1 ۱۵ ۰۲ ۱ 


نشريه مهندسی عمران فردوسی 


مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 
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شکل ۵ نمودار پراکندگی سرعت متوسط در عمق در اطراف 
خط ایده‌آل برای سه نوع شبکه انتخابی 


شرایط مرزی 

میدان جریان شامل توزیع سرعت و تنش برشی. 
باا ستفاده از مدل آ شفتگی ع- و نرم‌افزار -۸(5۷5 
و با دقت ۱*۱۰ حل شده‌است. شرایط مرزی 
به قرار زیر است: 

() توزیع یکنواخت سرعت در ورودی فلوم (۲) 
شرایط توزیع فشار هیدرواستاتیک در خروجی. (۳) از 
آنجایی که بدن؛ٌ فلوم از بتن صاف ساخته شده‌است؛ 
زبری دیواره‌ها و کف فلوم ۰/۱ میلی‌متر انتخاب 


شده‌است. 


سال سی و دوم شماره یک» ۱۳۹۸ 


محمدمهدی سیف- بهرام رضائی 


ارائة نتایج 

توزیع عرضی سرعت متوسط در عمق 
توزیع عرضی سر عت متوسط در عمق برای 
آزمایش‌های سری ۸14 و ۸15 بااستفاده از مدل 
آشضفتگی عع در دو مقطع ۱ و ۲ محاس‌به و با 
داده های آز مایش‌گاهی در همان مقطع مقایسه 
یله منت رفک هی( و (۷ توا پل 


د 


(ولمه) ۲۲۵ ,رهم۷۲۵1 
رت ه‌ِ 
مْ ح 
, 


۲ )۳( 
)2( 66۱100 1 


۷۵۱06, ۲۲۵ )/٩( 


۲ )۳0( 


)0( 66108 2 


شکل 7 توزیع عرضی سرعت در مقاطع ۱و ۲ برای کانال 


مرکب مورب آزمایش سری ۸14-2 


(و/۲۲۵)۵ رته‌ماه۷ 


(ظ) ۲ 


)2( 66۱100 1 


(٩/)و۲۷‏ رله۷۵10 


۲ )۳۲( 
)0( 60 2 


مرکب مورب آزمایش سری ۸15-2 


شالن سی و دوم شماره یک. ۱۳۹۸ 


.0 سس 


۱۷ 


نتایجی که از برر سی شکل‌ها می‌توان گرفت به 

قرار زیر است: 

۱. انطباق نسبتاً قابل قبولی بين توزیع سرعت اندازه- 
گیری‌شده و داده‌های به‌دست‌آمده از شبیه‌سازی 
عددی وحود دارد فقط درمورد آزمایش ۸15-2 
مدل آشفتگی ع- توزیع سرعت در سیلاب‌دشت 
همگرا را در مقطع شمارة ۱ زیاد. و در کانال اصلی 


به‌مقدار جزئی کم برآورد کرده‌است. 


وا کر ات ۱ رامق اوساات دش سس 
راست (همگرا) است که این امر می‌تواند ناشی از 
تغییر شکل هندسی کانال باشد که سبب خروج 
جریان از سیلات دستا همگرا و ورود آن به 


سیلاب‌دشت واگرا می‌شود. 


۴۳ سرعت‌های ماگزیمم به جای و سط کانال اصلی در 


۱ فصل مشترک کانال اصلی و سیلابدشت 
واگرا:(ذریافت کته جریان) اتفاق افعاده و تقریبا 
۸ تابور کت ان معوسستط سر تفر این 


8 با افرانشنی تامنه ارت مسلاف‌دشنت‌ها از 0/۱ جحه 


به ۹/۲ درجه اختلاف بین داده‌های آزماد شگاهی و 
شبیه‌سازی عددی افزایش می‌یابد. 
هم چنین مقادیر میانگین قدر مطلق خطای نقطه‌ای 


بین داده‌های آزمایشگاهی و نتایجم حاصل از شبیه‌سازی 
فلدی ااسفادم. از -رلنظه (ره۱) سعاسته ‏ شک‌است: 
مقدار متوسط این خطاها برای آزمایش‌های سری ۸14 
در مقاطع اندازه گیری ۱ و ۲ به‌ترتیب ۷/۶ درصد و 1/4 
درصد و برای آزمایش‌های سری ۸15 در همان مقاطع 


به ترتیب ۷/۵ و ۵ درصد می‌باشد. 


۱ 0100 


۳ »میم 9۵۷۷ 


نشریه مهندسی عمران فردوسی 


۱۸ 


در این رابطه هه ۷۷ سرعت به‌دست آمده از 
شبیه‌سازی عددی و وه ۷۷سرعت اندازه گیری‌شده و 


تعداد داده‌ها می‌باشد. 


توزیع تنش برشی مرزی 

تنش برشی مرزی در کانال‌ها اغلب در آنالیز سیستم- 
های رودخانه صرف‌نظر می‌شود. اما توزیع تنش برشی 
مرزی به‌دلیل این که در کالیبره کردن مدل‌های پیشرفته. 
در مطالعهٌ تعادل نیروها و نیز در محاسبات مربوط به 
انتقال رسوب کاربرد دارد. از اهمیت قابل ملاحظه‌ای 
برخوردار است. تنش برشی می‌تواند به‌عنوان معیاری 
برای تعیین محل سلول جریان ثانویه و نیز محل 
فرسایش و رسوب‌گذاری در کانال‌های طبیعی به‌کار 
زود 

بااستفاده از مدل آشفتگی ع- تنش برشی مرزی 
در کانال اصلی و سیلاب‌دشت‌ها برای دو زاوية اریب 
۱ و ۹/۲ درجه محاسبه و با مقادیر آزمایشگاهی 
مقایسه شده‌است. نتایج حاصل از این شبیه‌سازی برای 
آزمایش‌های سری ۸14-2 و ۸15-2 در شکل‌های (۸) 
قرف فان کاده قاس رین گنفن‌های: ترش 
مرزی نشان داده‌شده در این شکل‌ها دلالت بر وجود 
نقطهٌُ پیک در فصل مشترک بین کانال اصلی و سیلاب- 
دشت‌ها سمت چپ (سبلاب‌دشت واگرا) دارد. این 
نقطه ماگزیمم دارای مقداری به‌انداز؛ تقریباً ۲/۵ برابر 
تنش برشی متوسط در کانال است. شدت تنش برشی 
فا تسوت یگ کرادت بهسک | 
تقریباً ثابت است. مشابه اين پدیده در توزیع عرضی 
سرعت نیز مشاهده شده‌است. در مورد تنش برشی 
مرزی نیز با افزایش زاوية اریب از ۵/۱ درجه به ٩۹/۲‏ 
درجه اختلاف بین داده‌های آزمایشگاهی و نتایج 
حاصل از شبیه‌سازی عددی افزايش و دقت مدل 


آشفتگی 6-6 در شبیه‌سازی کاهش می‌پابد. 


نشريه مهندسی عمران فردوسی 


مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 
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شکل ۸ توزیع تنش برشی مرزی در مقاطع ۱ و ۲ برای کانال 


مرکب مورب آزمایش سری ۸14-2 


18 
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شکل ٩‏ توزیع تنش برشی مرزی در مقاطع ۱ و ۲ برای کانال 


مرکب مورب آزمایش سری ۸15-2 


مقادیر میانگین قدر مطلق طای نقطه‌ای بین 
داده‌های آزمادشکاهی و نتایج حا صل از شبیه سازی 
عددی تنش بر شی مرزی نیز باا ستفاده از رابطة (۱۱) 
محاسبه شده‌است. مقدار متوسط این خحطاها در مقاطع 


ا ان کی ۱ و ۲ جرا از ماش های نتب ۸۱۸ 


سال سی و دوم شماره یک» ۱۳۹۸ 


محمدمهدی سیف - بهرام رضائی 


به‌تردیب ۱۲/۱ درصد و ۱۰/۵ درصد و برای 
آز مایش های سری ۸15 در همان مقاطع به‌ترتیب 
۱ و ۱۶/۲ درصد می‌باشد. 


۰100 ۲6۵ صنا ۲ 
20 


> رات 0 
۱1( ۸ 170( 


در این رابطه سس« تنش برشی مرزی فقنیت الم 
از شبیهسازی عددی و وه تنش برشسی 
اندازه گیری‌شده و لا تعداد داده‌ها می‌باشد. 


نیروهای برشی ظاهری 

با مطالعة نیروی برشی ظاهری در فصل مشترک قائم 
بین کانال اصلی و سیلاب‌دشت‌ها می‌توان به مطالعةٌ 
اندرکنش بین کانال اصلی و سیلاب‌دشت‌ها پرداخت. 
به این منظور و برای تعیین نیروی برشی عمل کننده در 
هر ناحیه از کانال» محیط تر شده به تعدادی المان 
تفسیم و مطابق شکل ۲( نام گذاری شده است. تنش 
برشی محاسبه شده برای هر یک از المان‌ها به‌صورت 
عددی انتگرال گیری می‌شود تا نیروی برشی در واحد 
طول کانال به‌دست آید (رابطةٌ ۱۲ را ببینید). 

(۱۲ 7 12-1 , زنط رز ر5۳ 


که در این رابطه ٩۳‏ نیروی برشی در واحد طول 
ع مل کذ نده بر ناحية 1 ام 5 توزیع تنش برشی 
عمل‌کننده بر ناحیهٌ مورد نظر (1ام) و :۸ جزء طول 
در راستای محور عرضی کانال (محور ) می‌باشند. 
حال مایق شکل (0۰ در سیلابدقت همگرا 
و بین مقاطع ۱ و ۲ حجم کنترلی به‌طول م1 (۱۰ متر) 
در نظر گرفته و معادلة اندازة حرکت را به‌شسرح زیر 
برای آن تنظیم می‌شود. 
بک(ر۳٩‏ + م66۳۲ - ۷۷5 + و5 - بم۳ - بج۳ 
رصم - مورصمل + مور(صم) ع و(که)- 
۱۳( 
قن. یم وانظه. ‏ تسیز هام قفا 
هیدرواستاتیک وارد بر مقاطع شمارة ۱ و ۰۲ 5 مولفة 


سال سی و دوم» شماره یک. ۱۳۹۸ 


۱۹ 


نیروی عکس العمل دیواره درراستای محور 2 5۳6 و 
5۳7 متوسط نیروهای برشی در واحد طول کانال وارد 
بر کف و دیوارة سیلاب‌دشت سمت راست 
(سیلاب‌دشت همگرا)؛ تمرم و مع(صما۷) اندازه 
حرکت سیال ورودی به حجم کنترل و خروجی از 
حجم کنترل در امتداد محور 2 عع(۷]0) اندازه 
حرکت سیال خروجی از فصل مشترک بین سیلاب- 
دشت. و کانال اصلی در امغداد مرو و ت(قه) 
نیروی برشی ظاهری در فصل مشترک بین کانال اصلی 
و سیلاب‌دشت همگرا است. از آنجایی که عمق جریان 
در طول کانال اریب ثابت است. بنابراین مجموع نیروی 
فشار هیدرواستاتیک وارد بر مقاطع شمارة ۱ و ۲ با 
موف نیروی عکس‌العمل دیواره درراستای جریان برابر 


و در خلاف جهت یکدیگر هستنده درنتیجه مجموع 


این سه نیرو برابراست با: 


(۱ 0< ۲۲ - مو۲ - بو۳ 


بنابراین معادله پس از ساده‌سازی به صورت زیر 
درمی‌آید. 
< و(۸5۲) - بارر5۳ و6 و۷۷5 


و(0۷۵۵) ح بو (1۷۲۵۴) + وچ (ظ۳۵) 
۱( 


شکل ۱۰ نمای بالای حجم کنترل بین مقاطع شمارة ۱و ۲ در 
سیلاب‌دشت سمت راست (سیلاب‌دشت همگرا) 


نشریه مهندسی عمران فردوسی 


جملات مختلف رابطه (۱۵) از روابط زیر قابل 


محاسبه خواهند بود. 


۱31۱( 0 < ,۷۷۲5 
)۱۷( ۶ -< ب(۷/0۳) 
)۱۸ ۷ 2 و (۷۵۴) 
(۱۹( 02(۷۷ - ,۵)0 < ممو() 


با جای گذاری روابط فوق در م عاد لا (۱۵) و 
دوباره‌نویسی خواهیم داشت: 
+با(5۳ + ع5۳۲) - وکجعم < و(۸۹۳) 


(02 - ,)۵ - (و62021۷ - 0)8:0:1۷ 

)۲۰( 
در این رابطه. م جرم مخصوص سیال. ۵ شتاب 
۱۵ سح مان ی تحص ورن مس ست ٩‏ 
شیب طولی کانال:8 و 82 ضریب تصحیح اندازة 
حرکت سیال ورودی و خروجی. ,0 و 02 دبی ورودی 
به و حروجی از سیلاب‌دشت سمت راست. ,۷۷ و ۷۷ 
سرعت متوسط جریان در مقاطع شمارة ۱ و ۲ ۷۷۵۶۶ 
تفیل مرگ پین ان تما ی 
سیلاب‌دشت است. لازم به‌ذکر است که اندیس 3 
معرف پارامترهای مربوط به سیلاب‌دشت سمت راست 


اسبتیت. 


۱ 
۱ 
۱ 
1 
۱ 
سّ ۱ 
۱ 
۱ 
۱ 


(وه 3 د9) 1 
1۳ 


مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 


برای حجم کنترل مشابه در سیلاب‌دشت سمت 


چپ (سیلاب‌دشت واگرا) می‌توان نوشت: 


+با(و5۳ + ر8۳) - و08۷5 < ر(که) 


03(۶ - 2)04 + (,6۸0۸1۷ - و0)85091۷ 
,۲۱( 
در رابطهه فوق ۷۰ حجم سیال در حجم کنترل 
سمت چپ. 5۳ و و5۳ نیروی برشی متوسط در واحد 
طول کانال وارد بر دیوار و کف سیلاب‌دشت سمت 
چپ. :8 و 8۰ ضریب تصحیح انداز؛ حرکت سیال 
ورودی و خروجی. :0 و 04 دبی ورودی به و خروجی 
از سیلاب‌دشت سمت چپ. ۷۷ و ,۷۷ سرعت متوسط 
جریان ورودی به حجم کنترل و خروجی از حجم 
کنترل سمت چپ ۷ سرعت متوسط در فصل 
تشر کانین کافال افیلی ی اوقت اسه: 
باا ستفاده از معادلات (۲۰) و (۲۱) نیروی بر شی 
طامری وز قضسل مفسفرکه پین. کال میسن ز 
سیلاب‌دشت‌های همگرا و واگرای کانال مورب با 
زوایای اریب ۵/۱ درجه و ۹/۲ درجه محاسبه شده و 


در حدول‌های ۲( تا (۵) اووفه یله استت: 


شکل ۱ نمای بالای حجم کنترل بین مقاطع شمارة ۱ و ۲ در سیلاب‌دشت سمت چپ (سیلاب‌دشت واگرا) 


نشريه مهندسی عمران فردوسی 


سال سی و دوم» شمارة یک» ۱۳۹۸ 


محمدمهدی سیف- بهرام رضائی 


جدول ۲ جملات مختلف معادلهة اندازءٌ حرکت برای سیلاب‌دشت همگرا کانال مرکب اریب با زاوية انحراف 


)۸( 
)( 


--68 
۶ 
۱ 


-/۱ 


و۳ ضیالانی ماش ها هبتر زر کب اع سوت یی گرا کانان میک استب با واونه 


)۸۲( 
)( 


-۲ 
۳/۲۰ 
-7 


--۷۱ 


ول ۶ مراد مختافت: معادله اتقازه هرک یرای متلاب دشت هسگرا کانال مرکب ارشت. با رامیة 


(طه) 
(0) 


۵ 
۵ 
-۱۱/۳۹ 


-۹ 


دول ۵: ,مات متلف معادلة اندازه: حرکت بوای» بیلات فش واگرا کانال رکب اریب با-واویة 


)۸۵۲( 
)( 


-۹ 
۲۰ 
32۳ 


-۲ 


(۷) م(۷]۵8) 


۳/۳۵ 
۹/۱ 
(0 


۱۱۳/۹ 


(۱) م,:(۱۷۲۵۵) 


۳/۸۷ 
۱۳/۷ 
۰۷/۳ 


۱۳2 


() ءم(هم۳) 


9۹۹ 
۳۳/۹۸ 
۱۱۱/۵ 


(۱۱۲۱۱۸۱ 


() ءءبرصهم۳1) 


۷/۷ 
۳/۳۹۹ 
۱۳/۳۳ 


۳۱۳۱/۳۵ 


مه (ظ۷]۵) رز (جم]۷() 
() 


-:۳ 
۵ 
۹ 


-۲۸ 


)( 
۳/۳۹ 
۸324 

1/۹ 


(۱-1: 


۸(ظ۷۱۵) یز( جم]۷() 
() 


-۷۱ 
-۹ 
- ۰ 


۳ 


۸ (۱۷]۱۵) +( مطم]۱۷) - 
() 


1۱۷ 
۲/۹ 
۱۹/۷۲ 


۳۳۱/۹۹۵ 


با بررسی و مقايسة جدول‌های فوق می‌توان نتایج 

زیر را استنتاج کرد: 
۱ مقادیر منفی نیروی برشی ظاهری در فصل مشترک 
ی تال مین یی هایگ 
واگرا بیان‌کنندة این واقعیت است که جریان سریم‌تر 


شسالن سی و دوم شماره یک. ۱۳۹۸ 


م2 (ظ۷]۱۵) + ی( ما ۷) - 


(۵) بارر5۳ + م5۳۲) 


۸۳۸ 
۱/۳ 


۳۳/۳۳۹ 


۱۳۳۷/۱۸ 


بارو۳٩‏ + ر5۲) 
(۲) 


۱۰/۷/۸ 
۱۸/۸۵ 


۳۵/۰ 


29۰/۹ 


(0۷) با((5۳ + م5۳) 


۱:۶ 
۱۸/۹۲ 
۱۳: 


۰/۸۷۱ 


بلرو۹۳ + ر6) 
() 


۱۱۱/۳ 
۱/۳۷ 
۳۹/۷۸ 


۳/۸ 


۱۷۹ 
)( 


1/۱۳۱ 
۷/۹۹ 
۱۷/۵ 


۳۹/۳۰ 


۱ 
)( 


9/۷ 
۱/۳ 
۳۳/۲ 


۳۳/۰ 


۱ 
)( 


۹۲ 
۱۱۱ 
۳۳ 


۱۳-۱ 


۳۷۹ 
)( 


۳/۹۳ 
۱۷/۱ 
اور 


۳۳/۸۹۵ 


۱ درحه 
1 ۱۳۹ 
۰/۵ ۸14-1 
۰/۳ ۸142 
۰/۱ ۸14-3 
۱/0۰ ۸14-4 


نحراف ۵/۱ درجه 


۱۹ ۱۱/۹ 
۸14-1 ۰/۵ 
۸142 1234 
۸14-3 ۰/۱ 
۸14-4 ۰/0۰ 


نحراف ۹/۲ درجه 


۹ 12 
۸15-1 20 
۸15-2 ۰/۳ 
۸15-3 ۰/۱ 
۸15-4 ۰/0۰ 


نحراف ۹/۲ درجه 
1۳21 ۹ 
۰/۵ ۸15-1 
۰/۳ ۸15-2 
1/3 ۸15-3 
۰/0۰ ۸15-4 


۳۹ 


در کانال اصلی. جریان در سیلاب دشت‌ ها را شتاب 


1 


۲ با افزایش عمق جریان عموماً قدر مطلق مقدار 


نیروی برشی ظاهری در فصل مشترک بین کانال 


نشریه مهندسی عمران فردوسی 


۱۹۲ 


۳ با افزایش زاوية اریب سیلاب‌دشت. عموماً قدر 
مطلق نیروی برشی عمل‌کننده در فصل مشترک 
افزايش می‌یابد که به مفهوم افزایش اندرکنش 
جریان بین کانال اصلی و سیلاب‌دشت‌ها است. این 
پدیده می‌تواند ناشی افزایش تبادل جرم و به‌دنبال 
آن افزایش تبادل اندازٌ حرکت بین سیلاب‌دشت‌ها 
و کانال اصلی باشد. نتایج مشابهی در کانال‌های 
مرکب غیرمنشوری با سیلاب‌دشت‌های همگرا 
نقویستاه شبات ق نات [17] سخاهته قنل‌اس: 

حال می‌توان بااستفاده از رابطة زیر نیروهای 
برشی ظاهری را به‌صورت درصدی از نیروی برشی 
کل مقطع کانال محاسبه و تغییرات آن را برحسب عمق 

نسبی ترسیم نمود (شکل‌های (۱۲) و (۱۳) را ببینید). 


کم 


* 0 (۲۲ 
28۸: 


9۸۹۳ < 


در این رابطه ۸٩۲‏ نیروی برشی ظاهری: ,۸ 
مساحت مقطع عرضی کل کانال و ,5 شیب کف کانال 


نت 


‌ ۱ ف‌ 
۰ 7 ل 
1 
23 
1 ۰4 
| 
۲ 9۸۹۲۳۱ 
5 
ِ ( ۸۹ 
1 (و۸۹ ۱ 
و 
15 5 0۲ 53 15 25 


شکل ۱۲ مقايسهٌ درصد نیروی برشی ظاهری محاسبه‌شده 
پااستفاده از مدل ع-ع و داده‌های آزمایشگاهی در کانال مرکب با 


سیلاب‌شت‌های مورب ۵/۱ درجه 


نشریه مهندسی عمران فردوسی 


مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 


050 ۰ 
040 4 
22 0.30 
ِ 
« 
1 020 ۱ 
5 .0۸۹۲ 
010 ( م۵۸۹۲« 
۸ 
000 
15 5 5 15- 25 


۸۹ 
شکل ۱۳ مقایسة درصد نیروی برشی ظاهری محاسبه‌شده 
بااستفاده از مدل ع-ع و داده‌های آزمایشگاهی در کانال مرکب با 


سیلاب‌شت‌های مورب ۹/۲ درجه 


نتیجه گیری 

در این تحقیق میدان سرعت در یک کانال مرکب با 

سیلاب‌دشت های مورب با زوایای ۱ درحه و ۹/۲ 

درجه بااستفاده از مدل آشفتگی ع-» شبیه‌سازی شده 

و نتایج با داده‌های آزمایشگاهی مقایسه شده‌است. 

مهم‌ترین نتایج حاصل از این شبیه‌سازی به‌شرح زیر 

اتتاه 

گر ع-ع می‌تواند توزیع سرعت را با دقت 
قابل قبولی شبیه‌سازی کند. 

۲. به‌دلیل تغییر شکل هندسی سیلاب‌دشت‌های کانال 
مرکب اریب. سرعت حداکثر در نزدیکی فصل 
شش کت شیم کانان لیف ادا دیق | گرا 
اتفاق می‌افتد. 

۳ مقايسة توزیع تنش برشی مرزی شبیه‌سازی‌شده با 
مقادیر آزمایشگاهی برای کانال مرکب با سیلاب- 
دشت‌های مورب با زوایای اریب ۵/۱ درجه و ٩/۲‏ 
درجه نشان می‌دهد که انطباق قابل قبولی بین آنها 
برقرار است. اما شبیه‌سازی تنش برشی درمقایسه 

۶ نیروهای برشی ظاهری در فصل مشترک عمودی 
نییان ۵ انتت تا علامت منفی نیروهای برشی ظاهری 
بیانگر اثر افزایندگی جریان کانال اصلی برروی 


سال سی و دوم شماره یک» ۱۳۹۸ 


محم1ا مهدی سیف بهرام رضائی ۱۳۳۳ 


لیر رین ون اشیاا دش هقاس ناد فصل مشترک قائم بین کانال اصلی 


و 


دشت‌ها عموما قدر مطلق نیروی برشی ظاهری در افرایش اند رکتش,ننن این قسفت‌ها است. 


۳ 


0۱۵ ,۳۱۵۱۷۵ موجن ممامموصم صز ماگ اصملناهانها] .ریز رتا2ع۲ 220 نه بهع0۳0108 1 .. 


199(۰) ,21-41 .۲0 ,117 ۷۵۱۰ ,تآن)فظ ,و۵0 و۱ ۳۵۲۵۱۵۱6 


۹ ,۲۱۵۱۷۰ متفام۳۱۵۵0 ماما عقطا فعمامصومنوظ متتامجممی رل .ظ عمط مه ۷ رومصصهژ .. 


,۱۵۵106011108 0۲ وم زصنخ دنا رموصواهنز ۵۶ ]مممنهموو نا تمهت 0 طاقن متلنتهیل زا رجمتامع۹ ممصدتون(۲ 


1977 


باماکتظ ۶ه انوه۲نصنا م۲ رعتعمط] ساطظ ,۲۱۵۲۷۰ تمصصفت سمل ماصا ممتاهعتاوه1۳۲ مه ره ب) ٩.‏ باامتاا . 


1990(۰) 
0۵ ۳۵۵۵605 ۲۱0 91607۵0 :هام۳2 امصصجن ۳۱۵۵0 ۳۲و رل ۲۰ 5 رصتااهگ 20 بظ ‏ .9 با0تل۲1 
1990(۰) ,197-214 ۳0۰ ,28 ۷۵۱۰ مکاظا بواع ۳۲۵۵/۱6۵0 0 


مهاد-0۷0 آطونمند بخ مه ماه لقاممصتومنظ رچنلزمد۳ امصصع ۳۱۵۵0 بکلط۳و. ول ب۲ نک مطلللهگ . 


(1993) ,1 ,۷۵۱ ,ا0افترظ ۶ تیه تلا م1 ر فمتتم‌مهه صتقام۶۱000 ٩16۵/۵۵‏ رعامصصعطه 

2 20 9۳۱۷۲ عمط هصملونا ععمصمطمام مزونه ۷ طازه مامصصفق متامصوزظ ماورصته ه عصتاام۱۷۵0 ر.ل یلهاعلط) 
و۲666 تاوظ رعمتهام۳۲۱۵۵0 ماه از وامصصف ‏ متامصصوز مهملن ‏ طاممصی. ما ما کف ناه 
2009(۰) ,۲81۷۵۲۹1۲7 حصحطعصتصص را[ 

0 0۱۱۵ ۱۷۲۵0011882 [2مزتمصتا تودا ومز۴ مصصوت مومت ما ممنا موز کالهم۵۱ ۲۷ و.ت) ب8ع۳۵221۴8 
1994(۰) ,1176-1197 0۰ ,120 ۷۵۱۰ ,۱6۵178ع ع۵ 3 

مهم صد آ۲ممرکصه۳ ]۰ اممحصنلع قمع م۲۱ ۵۶ عمتاام۷۵0 لهمزتعصصتا ناو ر.] رقعلمه ۷۷ 220 ۷۷۰ رتم ,۷۷ ,۷۷۷ 
2000(۰) ,4-15 .۲0 ,126 ,۷۵۱ رفظ ,و۵۱ ۵ ول عامصصفت 

مهو ۵ مصتاام۱۷۱۵0 آهمزمصصتا۲ نلق. ربظ نک لا بهعا 4ص2 .۲ روعفقظ ر.] 9۵6950 ب. ۳۲1506۲-۵026 
-303 ۴0۰ ,39 ۷۵۱ مکگضا ,واع۵0۲ ع 0 ۵۴ ۵ طمتامامعع ۷ 0مومموطاناو ان سمل امصصق) 
(2001) ,310 


امصصحت مصتاممومصمت ۵۶ ممته‌انصتی لمعتمصتاط لمممتعصمصصزنا معط ]۰ ریق .ظ رصملعاظ م2 .ظ رصهنه ۲ طاومصصره؟ :10 


2003(۰) ,645-652 ۲۵۰ 129 ,۷۵۱ ,ناف بع۵ ۵۴۳۱۲۵۱ ۵ ول , وما۳۲ 


۵۲ 20 طاز زورما آمصمهطمممن0۵ 4ماممومصم ۵۲ عصتام۱۷۱۵0 مممملباطا نت 1 تا .9 رتم 220 ,۲۲ رعصهک :11 


۷۷۵۱۵۲ 17 ۸0۷۵۳۱۵۵۶ ,۱۷۲۵061 ووع90 108مطوم۴! معط عصتونا صنقام۲۱۵۵0 عمط مه ممتاماهعع ۱۷ 
.(2006) ,1650-1664 .0 ,29 


۶ زنل زان رولوعظ1۳ حاوظ ر"عامصصقی مه وامانا بهلدم‌تژممما ۶۵ وزم۱۷۲۵0 مممملتاطانت 1 .نا زامن :12 


] 0000, )1993( 


۹ ۲0۲ 0۲۲ 0۲ طمتا2112 ۷ 60انما1(۵ رآ رتعوعع9)0 24 .ظ ,۲۱ رطه۱۷]۵۲۲ رگ .۸۵ ,امامووم ریت ٩.‏ باطع ۷۷۲ :13 


سال سی و دوم شمارة یک ۱۳۹۸ نشریه مهندسی عمران فردوسی 


۱۹۶ مطالعةٌ عددی اثر زاویة اریب د رکانال م رکب با سیلاب‌دشت‌های... 


2004(۰) ,۲6217 ,عمام2 ,۲/۵۷۷ ۳۶۷۵ ,"5امصصعط ]طهنهای م1 

تلود رد ممتمصتاو مممه۵ وم ت۲0 دما امصصق مهو ۵۶ ممتاقاتصتوی ۲00 موتقا. ر.۲ رصهصصههظ :14 
2010(۰) ,۱۲۲۵۱۷۵۲51 حصفطع تاه لا روروعط 1" 

-۱0 ۳۳1 وامصصقت مصتامموصصمن صا ۳۱۵10 ما ۵۶ مصتام0م]۷ امعتنممت ۲ بب۲ رتتتصه 20 .ظ 5626 ,15 
.(2018) ,2413-2424 .0( ,25 ۷۵۱۰ م4 ۳۵۱۱۵6۵ 5۵۵۵۵۵ ]0۵ 0۱۱۵ , فصتهآ۳۱۵000 عتاقصوزتم 

0 مارا زاصته)تهعصنا ۴ هصتادممم؟ مج ممتاقصصتای؟ ز۵ مع۳۵6۵0ظ۳ و.ل .ظ رمطمجم۴ م2 .نا بقتطات ر.ظ ,1 بکلتامی .16 
2008(۰) ,1-4 .0 ,130 ۷۵۱ ,اف ۱60و 0۳۱۵ 0۵ , عجمای‌تانن۸ نا ما حمتلع 115016002 

عقوت صما طزس عامصصقت 0صمممصمیی ما ما هام۵ .۷۷۰ .نا ,کتک مه .ظ ۳۴622 .17 
2011(۰) ,815-824 ۲0۰ ,۷۵۱/۱37 ,یآتافظ رق0ع۱۱و ع۱ ۵ ]۵ آمسول , فصنهام ۳۲1۵00 


نشريه مهندسی عمران فردوسی سال سی و دوم شمار یک. ۱۳۹۸ 


